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摘 要 : 针对 现 有 算法 对 图 像 边 缘 细 节 增 强 不 足 及 无 法 有 效 控制 各 尺度 信息 增强 程度 的 问题 ， 提 出 了 多 级 分 解 的 
Retinex 低 照 度 图 像 增强 算法 。 该 算法 在 Retinex ee et 通过 设置 不 同 的 滤波 参 才 ， 获取 表 
征 图 像 不 同 尺 度 信息 的 反射 分 量 和 照度 分 量 ; 通过 使 用 指数 函数 对 分 解 得 到 的 各 级 反射 分 量 进行 增强 ， 能 够 有 效 提 
升 图 像 边缘 细节 的 表达 能 力 ; 通过 使 用 S 型 函数 对 最 终 的 RURAHEETAUS, uu d 图 像 整体 亮度 的 同 
时 抑制 高 亮度 区 域 ; 通过 颜色 恢复 函数 对 增强 图 像 进行 后 处 理 , 进一步 差 和 失真 的 问题 。 实 验 结 果 表 明 ， 
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Multi-level decomposition Retinex low-light image enhancement algorithm 


Wang Ping!, Sun Zhenming?! 
> (1. College of Economics & Management, Beijing Polytechnic, Beijing 100029, China; 2. College of Resources & Safety 
h Engineering, China University of Mining & Technology, Beijing 100083, China) 


Abstract: The existing algorithms can't enhance the details of image edges and effectively control the enhancement degree 
of each scale information. Aiming at the problem, this paper proposed a multi-level decomposition Retinex low-light image 
enhancement algorithm. On the basis of Retinex decomposition model and bilateral filtering, the algorithm can obtain the 
multi-level reflection components and illumination component by setting different filtering parameters. Using the 
exponential function to enhance the multi-level reflection components, the algorithm can effectively enhance the expression 
ability of the image edge details. And using the S-type function to process the final luminance component, the algorithm can 
improve the overall brightness of the low-light image while suppressing the high-luminance area. By processing the color 
restoration function as post-processing, the algorithm can further avoid the problems of color deviation and distortion. The 
experimental results show that the proposed algorithm can improve the visual quality of low light images and is superior to 
the existing algorithms in terms of clarity, information entropy, and contrast. 

Key words: bilateral filter; multi-level decomposition; low light image enhancement; Retinex algorithm 


0 引言 终 增强 的 图 像 。 对 于 低 照 度 图 像 而 言 ， 早 期 的 Retinex 相关 
x: 算法 存在 严重 的 色彩 失真 现象 。 为 解决 上 述 问题 ， 现 有 的 
随 着 信息 化 手段 的 发 展 ， 生 产生 活 中 产生 了 大 量 低 照 度  Retinex 算法 同时 考虑 照度 分 量 和 反射 分 量 对 图 像 的 影响 , 通 
境 下 获取 的 图 像 。 这 些 图 像 具 有 整体 亮度 低 、 分 辨 率 低 、 过 使 用 不 同 的 增强 函数 对 分 解 得 到 的 照度 分 量 和 反射 分 量 进 
噪声 大 等 特点 ,制约 了 人 眼 辨识 和 进一步 的 图 像 处 理 与 分 析 。 ” 行 处 理 从 而 有 效 提升 图 像 的 质量 ， 典 型 的 算法 包括 非 均 匀 照 
为 了 提升 图 像 可 视 化 效果 ， 有 必要 设计 专门 的 低 照度 图 像 增 。 度 图 像 色 彩 恢 复 算 法 隔 、 低 照度 图 像 增强 算 法 外 等 。 然 而 ， 
强 算 法 。 传 统 的 低 照度 图 像 增强 算法 可 以 分 为 基于 伪 雾 图 的 上 述 算法 在 处 理 图 像 时 缺少 对 图 像 多 尺度 信息 的 考虑 ， 使 得 
曾 强 算 法 、 对 比 度 增强 相关 算法 和 Retinex 相关 算法 。 由 于 ”算法 那 对 边缘 细节 等 高 频 信息 增强 不 足 ， 限 制 了 图 像 视觉 质 
轨 度 图 像 的 负片 与 雾 者 图 像 非常 相似 ， 基 于 伪 雾 图 的 增强 — 量 的 提升 。 
法 能 够 根据 雾 才 成 像 机 制 对 低 照 度 图 像 的 负片 进行 “ 雾 霆 ” 为 了 有 效 解 决 现 有 算法 对 图 像 多 尺度 信息 考虑 不 足 的 i 
晰 化 处 理 ， 最 后 转换 得 到 增强 后 的 图 像 上 23。 对 比 度 增 强 算 — 题 ， 本 文 提 出 多 级 分 解 的 Retinex 低 照度 图 像 增强 算法 。 该 
通过 调整 图 像 的 像素 或 直方 图 分 布 来 提升 图 像 整 体 的 对 比 ”算法 通过 迭代 使 用 双边 滤波 实现 对 图 像 的 多 级 分 解 ， 获 取 表 
， 实 现 低 照 度 图 像 的 增强 ， 典 型 的 算法 包括 受 限 对 比 度 自 。”” 征 图 像 不 同 尺度 细节 信息 的 照度 分 量 和 反射 分 量 ， 改 进 图 像 
立 直 方 图 均衡 化 算法 中 、 广 义 非 锐 化 掩 膜 算法 外 等 。 增强 函数 ， 使 用 对 数 增强 函数 和 S 型 增强 函数 对 获取 得 到 的 
Retinex 相关 算法 认为 图 像 可 以 表示 为 照度 分 量 和 反射 ” 多 级 反射 分 量 和 照度 分 量 进行 增强 ， 最 后 ， 融 合 各 分 量 ， 获 
量 的 卷 积 ， 其 中 照度 分 量 反映 整体 变化 缓慢 的 光照 信息 ， 得 最 终 的 增强 图 像 。 
反射 分 量 才 是 图 像 的 本 证 属性 。 早 期 的 Retinex 相关 算法 
通过 去 除 照度 分 量 对 图 像 的 影响 ， 将 获取 的 反射 分 量 作为 最 
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1.1 

它 能 
Retinex 到 
Retinex 算法 、 


的 Retinex 算法 、 基 于 偏 微 分 方程 的 阐 值 Retinex 算法 以 及 基 


Retinex 理论 及 存在 问题 
Retinex 图 像 分 解 模型 


滤波 与 Retinex 


于 变 分 的 Retinex 算法 。 


um 
结合 , 通 


1986 年 ，Jobson 等 人 四 将 高 斯 低 通 
过 改进 Land 等 人 外 提出 的 中 心 环绕 


Go 萍 ， 等 : 多 级 分 解 的 Retinex 低 照 度 图 像 增强 算法 


是 基 ZTA 民 视 觉 系 统 提 出 的 图 像 增强 理论 ， 
彩信 息 的 同时 有 效 增强 图 像 细 节 。 围 绕 
里 论 产生 了 一 系列 图 像 增强 算法 
基于 随机 采样 的 Retinex 算法 、 基 于 中 心 环绕 


,包括 基于 路 径 的 


Retinex 算法 (center/surround Retinex)， 提 出 了 单 尺 度 


Retinex (SSR) 算法 。 


Jobson 等 人 使 上 


不 是 由 反 
HJER 


可 见 光 图 像 的 场景 光照 信息 ， 获 取 反 射 
Xt ER EE 


mE ( 绿 ) 和 短波 〈 蓝 


ENS K, 


分 量 


i 低 通 滤波 获取 照度 分 
量 的 方式 使 算法 的 执行 效率 有 了 极 大 的 提升 ， 称 这 种 方式 为 
基于 分 解 的 Retinex 算法 。 
里 论 的 基本 内 容 是 物体 的 颜色 是 由 物体 对 长 波 
) 光线 的 反射 能 力 决 定 的 ， 而 
强度 的 绝对 值 决定 的 ; 物体 的 色彩 不 受 光 照 非 


致 性 ， 即 Retinex 理论 是 以 色 感 一 臻 
性 (颜色 恒 党 T Heb 为 基础 的 。 Retinex 算法 的 核心 是 通过 去 除 


< 图 像 的 反射 分 量 ， 
: 原 图 像 ，@ 表示 卷 积 操作 。 
通常 ， 物 体 的 反射 分 量 在 全 部 反射 和 全 部 吸收 之 间 ， 即 


反射 分 量 的 取 值 范围 满足 FG,y)e[0,1] 。 


分 量 满足 如 下 关系 : 


L(x,y)=U,y)/ FG, y) ZUG y) (2) 


1.2 Retinex 图 像 增 强 模型 


习 像 分 解 为 照度 分 


像 分 量 ， 从 而 实现 


型 ， 可 以 将 图 像 看 做 照度 分 量 和 反射 分 


Ux, y) FG, y) 8 Læ, y) (1) 
LE, y) 表示 图 像 的 照 


Kd 


进而 ， 可 以 推断 照度 


量 和 反射 分 量 , 分 别 刻 


H f B 


一 一 百 


r 


。 为 了 使 图 像 获得 如 


使 用 不 


a 


U‘, y) -f 


度 分 量 和 反射 分 量 进行 增强 ， 进 
BORSE NEN AE 和 反射 分 量 ， 获 得 最 终 增 强 后 的 图 
车 于 分 解 的 Retinex 图 像 增强 模型 


更 好 的 可 视 化 效果 ， 


ChinaXiv 合 作 期 刊 


E 


了 许多 文章 在 伽 马 函数 和 $ 型 函 
的 增强 函数 。 


WB Retinex 图 像 增强 流程 
Fig.1 Flowchart of Retinex based im: 


输入 图 像 
| 图 像 分 解 


各 图 像 分 量 增强 
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age enhancement process 


Kf] Retinex 算法 使 用 丰 


的 关系 。 为 此 ， 


度 Retinex 算法 


尺度 处 理 策 咯 的 关键 是 全 
使 得 获取 的 照度 分 量 包含 不 同 尺 


` 同 的 滤波 函数 来 获取 图 
于 滤波 尺度 不 好 把 \ 能 平衡 图 


映 图 像 不 同 尺度 信 
Æ 9 如 权 和 的 方法 不 能 较 好 反映 体 图 像 的 特点 ， 


使 用 不 同 


度 增强 程度 的 有 效 控制 ， 
边 滤波 的 基础 上 提出 多 级 分 解 的 滤波 策 


量 和 最 终 照 度 分 量 进行 增强 ， 实 现 对 
， 从 而 能 和 图 像 整 


4 可 以 表示 为 
Ea, y) 8 g LG y) 6) 


数 ，gO 表示 照度 分 量 的 增强 函数 。 基 
on 如 非 均匀 照度 图 
度 图 像 增强 、 细 节 增 强 、 


HR RZ, fO 表示 反射 分 量 的 增强 函 


模型 已 经 广泛 
像 增强 、 低 照 


于 分 解 的 Retinex 图 像 


高 动态 范围 图 像 映射 


等 。 针 对 不 同 包 


k 理 领域 前 需要 


针对 性 的 设计 相应 的 


2 DASS 理 需 求 。 型 


pa 


奋进 行 考察 ， 即 获取 照度 分 量 


EL 


的 增强 函数 以 及 多 尺度 信息 的 保持 。 


&| 1 所 示 。 


的 基于 分 解 的 Retinex 


基于 分 解 的 Retinex 算法 的 核心 要 素 可 以 从 以 下 三 个 方 
的 滤波 器 、 照 度 分 量 和 反射 分 


照度 分 量 反映 的 是 图 


像 变化 缓慢 的 信息 ， 


通 滤波 等 滤波 方法 。 


早期 研究 者 主要 使 用 高 斯 滤波 、 频 域 低 


在 模糊 pu 


。 近 几 年 


bh 


度 分 量 进 行 调 整 。 对 于 


实现 图 像 照 度 的 均衡 ， 即 对 低 照度 


于 这 些 滤波 方法 不 具有 保 边 特性 ， 即 
保持 图 像 的 边缘 ， 从 而 使 得 最 终 获取 
9 5、 迭代 中 值 滤波 、 
SERIAS AE REREN EI E 
于 基于 分 解 的 Retinex 算法 中 。 


， 双 边 滤波 上 


等 高 频 信 息 。 
多 级 分 解 的 Retinex 图 像 增强 模型 


虑 波 等 方法 导致 的 颜色 失真 和 光 晕 效应 ， 
滤波 、 引 导 滤 波 等 具有 边缘 保持 4 


尺度 信息 的 获取 ; 同时 通过 使 用 增强 函数 对 不 同 尺度 


的 反射 


VERRE 


息 的 照度 分 量 
尺度 照度 
同时 不 
为 此 本 文 在 双 
各 ， 实 现 对 图 像 不 同 


尺度 信息 有 效 控 


不 同 


的 
1$] Retinex 相关 算法 中 。 双 边 滤波 是 在 高 
能 


2n d E 分 析 
滤波 方法 在 图 像 处 理 中 有 广泛 的 应 | 


方面 调节 照度 图 像 的 生 


了 保 边 特性 的 滤波 方法 ， 同 时 


亮度 的 同时 增 


REPER 


CEH 


加 权 均 值 来 实现 图 像 光 ; 


“对 于 全 时 度 图 像 处 理 ， 可 以 使 用 伽 马 函 数 来 对 图 像 的 昭 


F 非 均匀 照度 图 像 ， 可 以 使 用 S 型 函数 
区 域 进 行 增强 、 对 高 照度 


J, CEARR 
并 减少 图 像 中 的 噪声 。 在 滤波 过 


W, -» G,, (lp m 4| )G,, (U, e U,) 
qe 


斯 滤波 核 的 值 只 与 空间 领域 的 距离 有 关 ， 而 与 图 像 
的 像素 值 无 关 ， 即 每 个 像素 点 的 滤波 核 都 一 致 ， 而 不 能 
的 图 像 上 下 文 进 行 调 整 ， 从 而 使 得 
和 性 。 为 了 解决 上 述 问题 
网 觉 大 会 上 提出 了 
。 双 边 滤波 的 公 


1 
Ul, = 2,6. (Ip -aqDGs, U, -UU, 
p qe 


根据 


时 获 取 的 模糊 图 像 
题 ，Tomasi 和 Manduch 
边缘 保持 特 
式 可 以 表示 为 


(4) 


(5) 
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其 中 : W 为 归 一 化 因子 , 使 图 像 的 像素 值 在 滤波 前 后 处 在 统 。 特征 ， 使 用 具有 较 大 滤波 半径 Cmn ) 的 双边 滤波 对 一 级 
一 范围 ，P 和 4 表示 图 像 中 像素 的 位 置 ， 为 当前 的 滤波 位 。 照度 分 量 工 进行 滤波 ， 获取 图 像 的 二 级 照度 分 量 工 ,根据 式 
置 ，? 表示 ? 的 令 域 9 中 的 像素 位 置 ，V, 和 0 分别 表示 像 ” (2) 对 照度 分 量 进行 修正 ,使 其 满足 懂 >L ; edis P = L/L 3R 
R P FI 处 的 像素 值 。 其 中 6。 和 Go, 表示 二 维 高 斯 滤波 核 ， 取 图 像 的 二 级 反射 分 量 F 。 


即 c) 按照 同样 的 方式 获取 更 高 级 的 反射 分 量 和 照度 分 量 ， 
Ga- agp) (6) Mes pe 

^5 229 c? 图 像 多 级 分 解 能 够 有 效 提取 图 像 不 同 尺度 的 信息 。 在 使 

XX (6) 可 以 看 出 ， 双 变 滤波 的 核 由 两 部 分 组 成 ， 即 用 的 过 程 中 ， 需 要 注意 如 下 几 点 ， 为 提取 图 像 不 同 尺 度 下 的 


G., (p -4 FI Ge, (U, -U,), Go, 由 空间 领域 距离 构成 ，G- 由 图 言 息 ， 对 原 图 像 滤波 使 用 的 滤波 半径 应 该 比较 小 ， 然 后 在 后 
像 的 次 度 差 值 构 成 。 相 应 的 oc 和 o 分 别 表示 相应 的 标准 差 ， 续 滤 波 中 持续 的 增加 ， 即 5<5<s<s<…。 不 同 于 传统 的 多 
一 般 而 言 ， 标 准 差 越 小 ， 平 滑 效应 越 不 明显 ;标准 差 越 大 ， 尺度 滤波 使 用 不 同 尺 度 的 滤波 器 对 同一 由 图像 分 别 进行 滤波 ， 
平滑 效应 越 明显 。 然后 取 滤 波 后 的 图 像 取 加 权 和 作为 最 终 输出 ， 本 文 提出 的 多 
于 双边 滤波 的 核 包 含 了 图 像 邻 域 像素 值 的 差异 ， 因 此 ”尺度 滤波 能 够 以 多 级 分 解 的 方式 获取 图 像 多 个 尺度 的 信息 。 

它 能 够 根据 图 像 邻 域 像素 值 的 差异 的 改进 而 自 使 用 调整 ， 进而 ， 图 像 可 以 通过 以 连 乘 的 方式 从 获取 的 各 尺度 反射 分 量 
实现 边缘 保持 的 目的 。 和 最 终 的 照度 分 量 来 重建 ， 这 为 控制 各 尺度 的 增强 结果 从 而 
图 2 展示 了 双边 滤波 不 同 参数 下 的 滤波 效果 。 适应 更 多 的 图 像 应 用 提供 了 可 能 。 通 过 多 次 实验 ， 发 现 使 用 

三 级 的 多 级 分 解 的 滤波 过 程 就 可 以 取得 较 好 的 效果 。 从 图 3 


可 以 看 出 ， 随 着 分 解 层 次 的 提升 ， 图 像 的 照度 分 量 越 来 越 光 
滑 ， 同 时 能 够 有 效 保持 图 像 中 的 边缘 ， 同时， 反射 分 量 所 能 
反映 的 高 频 信息 的 尺度 也 越 来 越 大 ， 即 反射 分 量 中 图 块 的 聚 
集 性 逐渐 提高 。 
| 
| 
| 
| 
| 
| 


图 2 ”双边 滤波 不 同 参数 对 图 像 模糊 的 影响 


(al) AR Cb1) 增强 图 像 


| 

| 

| 

| 

| 

| 

Fig.2 Influence of different parameters of bilateral filtering on image | 
| 

blur | 

图 2 中 可 以 看 出 ， 双 边 滤波 能 够 有 效 根据 不 同 的 滤波 | 
参数 实现 对 图 像 不 同 程度 的 模糊 ， 在 模糊 的 同时 能 够 有 效 保 | 
持 图 像 的 边缘 。 而 且 ， 随 着 参数 的 增加 ， 双 边 滤 波 对 图 像 的 | 
| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 

| 
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(a2) 一 级 照度 分 量 (b20 一 级 反射 分 量 


模糊 呈现 逐 级 递增 的 效应 ， 或 者 说 可 以 通过 使 用 不 同 参数 的 
双边 滤波 来 获取 模糊 程度 不 同 的 图 像 序列 。 一 般 而 言 ， 当 使 
用 小 尺度 的 滤波 函数 对 图 像 处 理 时 ,图 像 的 细节 保留 比较 多 ， 
结合 Retinex 模型 可 以 获取 小 尺度 上 的 反射 分 量 ， 进 而 可 以 
实现 逐 级 增加 滤波 参数 ， 获 取 不 同 尺度 上 的 模糊 图 像 作为 照 
EDE; 并 逐 级 获取 反映 不 同 尺度 情况 的 反射 分 量 。 受 此 启 
发 ， 本 文 提出 了 图 像 的 多 级 分 解 ， 并 基于 Retinex 原理 获得 
不 同 尺度 下 的 照度 分 量 和 反射 分 量 。 


(a3) 二 级 照度 分 量 (b3) 二 级 反射 分 量 


三 级 反射 分 量 
22 基于 双边 滤波 的 图 像 多 级 分 解 ee ns 
根据 Retinex 理论 ， 使 用 带 有 较 小 滤波 半径 的 双边 滤波 图 3 基于 双边 滤波 的 图 像 多 级 分 解 示意 图 
对 图 像 处 理 后 仍然 保留 了 大 量 的 细节 特征 ， 这 为 图 像 的 多 级 Fig.3 Multilevel image decomposition based on bilateral filtering 
分 解 提供 可 能 。 在 此 基础 上 ， 本 文 提 出 的 图 像 多 级 分 解 的 示 图 4 展示 了 图 3 中 第 一 列 各 图 像 中 第 100 行 的 提取 线 示 
意图 如 图 3 所 示 。 意图 ， 其 中 蓝 色 表示 原 图 像 像素 分 布 ， 红 色 表 示 一 级 照度 分 
本 文 提出 的 图 像 多 级 分 解 的 步骤 为 : 量 的 光滑 情况 ， 由 图 中 可 以 看 出 ， 一 级 照度 分 量 对 图 像 的 光 
a) 使 用 具有 较 小 滤波 半径 i 的 双边 滤波 对 输入 图 像 进 滑 效 应 较 弱 , 对 应 于 图 3(a2) 包 含 了 大 量 的 纹理 细节 ; 绿色 表 
行 处 理 ， 获 取 图 像 的 一 级 照度 分 量 工 ， 根 据 式 (2) 对 照度 分 量 。 示 二 级 照度 分 量 ， 是 根据 Retinex 模型 基于 一 级 照度 分 量 获 
进行 修正 ， 使 其 满足 五 > 五 ; WEAR E =L /WL 获取 图 像 也 取 的 ; 青色 表示 三 级 照度 分 量 的 提取 线 的 像素 情况 ， 同 样 根 
一 级 反射 分 量 五 。 据 Retinex 模型 基于 二 级 照度 分 量 获取 ， 其 对 图 像 的 光滑 效 


b) 考虑 步骤 a) 获 取 的 一 级 照度 分 量 工 包含 丰富 的 纹理 。 应 比较 强 ， 反 映 了 图 像 较 大 尺度 上 的 变化 趋势 。 
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一 # 一 原 图 像 


一 人 一 一 级 照度 分 量 
04| 二 级 照度 分 量 
x 三 级 照度 分 量 


0 100 200 300 400 500 


图 4 图 像 多 级 分 解 的 效果 示意 图 〈 以 图 3 为 例 ) 
Fig.4 Schematic diagram of image multilevel decomposition effect 
(Fig. 3) 
2.8 多 级 分 解 的 Retinex 增强 模型 
图 像 的 多 级 分 解 能 够 将 图 像 表 达 为 一 系列 的 照度 分 量 和 
反射 分 量 ， 该 过 程 可 以 表示 为 
U=F 9L 
=F, 9F, 9L, (7) 
-F,GF,G..GF,GL, 
Hp: n 表示 最 终 分 解 的 级 数 , Fus k e Un] RR k 205 2) 365, 
L keln KR k 级 照度 分 量 。 本 文 将 式 〈7) 表示 的 图 像 模 
型 成 为 广义 化 的 Retinex 模型 。 在 此 基础 上 ， 可 以 使 用 不 同 
的 增强 函数 分 别 对 各 级 反射 分 量 和 照度 分 量 进行 增强 ， 然 后 
将 增强 后 的 各 分 量 进行 处 理 ， 得 到 最 终 增强 的 图 像 ， 即 
U" =F) 8£(F,) 8- ORF) 8 gL) (8) 
其 中 :tf 表示 各 级 反射 分 量 的 增强 函数 ; B) 表示 n 阶 照 度 
分 量 的 增强 函数 。 针 对 不 同 的 增强 任务 ， 可 以 采用 不 同 的 增 
强 函 数 对 各 分 量 进行 处 理 , 针 对 本 文 低 照 度 图 像 增强 的 应 用 
本 文 将 分 别 使 用 指数 函数 和 S 型 函数 实现 对 各 级 反射 分 量 和 
照度 分 量 的 增强 ， 有 具体 公式 将 在 后 文 进行 叙述 。 
多 级 分 解 的 Retinex 增强 模型 的 推导 过 程 如 下 : 首先 ， 
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进行 处 理 ， 然 后 对 亮度 通道 进行 多 级 分 解 ， 获 取 表 征 图 像 不 
同 尺度 的 反射 分 量 和 照度 分 分 量 ， 为 提升 图 像 不 同 尺度 细节 
言 息 的 表达 能 力 ， 使 用 指数 函数 对 获取 的 各 级 反射 分 量 进行 
处 理 ， 为 提升 图 像 整体 亮度 的 同时 抑制 高 亮度 区 域 ， 使 用 S 
型 函数 对 获取 的 最 终 照 度 分 量 进行 处 理 ， 融 合 增强 后 的 各 分 
量 得 到 增强 后 的 亮度 通道 并 转换 回 RGB 颜色 空间 ， 对 其 进 
行 色彩 恢复 并 输出 。 算 法 流程 如 图 S 所 示 。 


低 照 度 图 像 


RGB->HSI 


2 级 反射 分 量 F， 
3 级 反射 分 量 F， 


指数 函数 


32 UR EA EL, 


HSI-»RGB 


融合 图 像 


色彩 恢复 


输出 图 像 图 像 


图 5 基于 图 像 多 级 分 解 的 低 照度 图 像 增强 算法 


Fig.5 Low illumination image enhancement algorithm based on 


image multilevel decomposition 
3.1 颜色 空间 转换 
为 提高 图 像 处 理 的 时 间 效 率 同时 较 好 的 保持 图 像 的 色彩 


根据 Retinex 原理 ， 对 原 图 像 进行 分 解 ， 即 UU = eL, IHE 
一 级 反射 分 量 和 照度 分 量 。 考 虑 到 一 级 照度 分 两 种 仍然 包含 
了 大 量 的 纹理 细节 ， 从 而 为 图 像 的 进一步 分 解 提供 可 能 。 为 
此 ， 对 一 级 照度 分 量 进行 处 理 ， 即 L =F 8L 获得 二 级 反射 
分 量 和 照度 分 量 。 依 此 类 推 ， 可 以 获得 式 CIO 的 表述 。 在 获 
得 各 级 反射 分 量 及 照度 分 量 的 表示 后 , 根据 式 (8) 对 图 像 进 
行 增强 。 其 中 需要 解决 问题 ， 即 增强 后 的 像素 值 范 围 不 能 超 


Tr 


ci 


cl 


EKE, 根据 Retinex 原理 ， 照 度 分 量 L(x,y) 反映 了 入 射 
到 物体 上 的 光源 能 量 总 和 ， 它 表示 了 图 像 中 变化 缓慢 的 低频 
Hi, 满足 LOS») el0,%] ,反射 分 量 F(,y) 是 物体 的 反射 系数 ， 
它 包含 了 图 像 的 高 频 的 细节 部 分 ， 满 足 7(x,y) el0,1] 。 通 常 ， 
物体 的 反射 分 量 在 全 部 反射 和 全 部 吸收 之 间 ， 即 反射 分 量 的 
取 值 范围 满足 F(x,y) el0,1] 。 
在 进行 图 像 的 级 联 分 解 过 程 中 ， 本 文 将 图 像 从 [0,255] 调 
整 到 [0.0,1.0] 进 行 处 理 。 这 样 处 理 的 好 处 是 能 够 最 大 限度 的 保 
证 增强 后 的 图 像 像 素 值 限定 在 [0.0,1.0] 内 , 从 而 不 会 造成 图 像 
的 过 增强 或 欠 增 强 问 题 。 事 实 上 ， 根 据 式 (7)， 本 文 获取 的 
各 级 反射 分 量 和 照度 分 量 的 值 都 限定 在 [0.0,1.0] 内 。 


3 ” 低 照 度 图 像 增 强 算 


为 了 对 低 照 度 图 像 进行 处 理 , 本 文 在 多 级 分 解 的 Retinex 
增强 模型 的 基础 上 提出 了 本 文 的 增强 算法 。 该 算法 首先 将 图 
像 从 RGB 颜色 空间 转换 到 HSI 色彩 空间 ， 并 提取 亮度 通道 


信息 ， 本 文选 择 对 HSI 颜色 空间 的 亮度 通道 进行 处 理 。RGB 
颜色 空间 转换 到 HSI 颜色 空间 可 以 表示 为 
I=(R+G+B)/ x43 (9) 
2R-G-B 
H -[90— f 
PECES agp) (10) 
+{0,G > B:180,G < B]]/360 
s (R-Gy «-(R- BYG - B) (11) 


对 于 HSI 空间 ， 各 颜色 分 量 的 取 值 分 别 为 0 <H «360 , 
1<S<1，1<1<1。 通 常 HSI 颜色 空间 到 RGB 颜色 空间 的 通 
常 按 照 色调 HH 的 不 同 分 区 进行 ,对 于 RG 扇 区 ( 0 € H «120 ), 
转换 公式 为 


B-I0-8) (12) 
ScosH 

= 60 8 (13) 

G —3I - (R+B) (14) 


对 于 GB AX (120 x H «246 ), $ H'=H-120 ， 则 转换 
公式 为 


R-I(1I-S) (15) 

S cos H! 
G-HI ev H (16) 
B 231 -(R«G) (17) 


对 于 BR HX 〈240 xH «3600, 4 H?-H-24* ， 则 转 
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换 公 式 为 
G -I(0-S) (18) 
B-diis S cos H? ] (19) 


cos(60° — H?) 
R-3I - (B +G) (20) 
3.2 图 像 多 级 分 解 和 增强 
在 获取 图 像 的 亮度 通道 之后， 使 用 本 文 提出 的 多 级 分 


解 对 其 进行 处 理 。 对 于 低 照 度 图 像 ， 首 先 对 其 进行 分 解 ， 分 
别 获取 表征 图 像 精细 的 小 尺度 信息 、 中 尺度 信息 和 大 尺度 信 
息 的 三 级 反射 分 量 和 最 终 的 照度 分 量 。 对 于 一 级 分 解 ， 设 定 x 
双边 滤波 参数 为 =5 ，o,=5.0 ，o,=5.0 ; 对 于 二 级 分 解 ， 设 图 6 不 同 参数 下 的 S 型 曲线 函数 
定 双边 滤波 的 参数 为 i=10 ，o, =100 ，o,=10.0 ; 对 于 三 级 分 Fig.6 S-type function under different parameters 
解 , 设 定 双 边 滤波 的 参数 为 1=30 ，o, 2150 ，o, 2150 。 最 终 ， 在 对 各 分 量 进行 增强 之 后 得 到 增强 的 亮度 通道 图 ， 即 
得 到 一 级 反射 分 量 瑟 、 二 级 反射 分 量 F,、 三 级 反射 分 量 包 和 Ie =f(F,F,,F,)® glL,) 24) 
三 级 照度 分 量 L. 然后 根据 式 (120 ~ QOO 将 图 像 由 HSI 颜色 空间 转换 
为 了 提升 各 级 反射 分 量 的 细节 表达 能 力 ， 采 用 指数 函数 可 RGB 颜色 空间 。 
对 各 级 反射 分 量 进行 处 理 。 即 3.3 颜色 恢复 
KF,.x,) = (F,)* Q1) 由 于 上 述 增强 方法 只 对 图 像 的 亮度 图 进行 处 理 ， 而 缺少 
其 中 : 下 标 k 表 示 反 射 分 量 的 序号 且 ke{l,2,3} ; 表示 指数 ”对 色调 信息 的 处 理 ， 容 易 造成 图 像 存 在 色彩 偏差 ， 为 此 参考 


增强 系数 。 由 于 具体 的 图 像 增 强 场 景 对 各 尺度 细节 增强 强度 。 带 色 彩 恢复 的 多 尺度 视网膜 增强 算法 (MSRCR) 对 图 像 进行 
的 要 求 不 一 样 ， 所 以 对 于 各 级 的 反射 分 量 应 设置 不 同 的 指数 MEKE. KRENE, HG G-12.3» 表示 原始 输入 的 低 照 
系数 。 本 文 提 出 了 基于 双边 滤波 的 图 像 多 级 分 解 ， 并 在 此 基 ERRAI R,G,B 颜色 通道 ， g; 0-12. 表示 待 色彩 恢复 的 图 
础 上 提出 了 多 级 分 解 的 低 照度 图 像 增 强 算 法 。 该 算法 的 核 ， f. E G-12.3) 表示 色调 恢复 后 的 输出 图 像 。 则 颜色 恢复 的 
是 通过 图 像 的 级 联 分 解 ， 获 取 表 征 图 像 不 同 尺 度 的 细节 信息 。 公式 可 以 表示 为 
的 反射 分 量 ， 进 而 在 对 使 用 不 同 的 增强 系数 对 各 级 反射 分 量 


E> 


进行 增强 。 u C, - f*ln 1e, 7G, Q5) 
在 2.2 节 ， 本 文 详 细 分 析 了 多 级 分 解 的 步骤 及 获取 的 各 G10+1 

级 反射 分 量 情况 。 一 般 而 言 ， 由 于 使 用 较 小 滤波 半径 的 双边 H - 

滤波 对 原始 图 像 进行 处 理 ， 其 获取 的 一 级 反射 分 量 表征 的 更 8, =C, 8g, (26) 

是 图 像 高 频 的 细节 信息 ; 随 着 滤波 半径 的 增 大 ， 双 边 滤波 其 中 : 8 表示 增益 参数 ; 7 用 来 调节 彩色 图 像 的 亮度 程度 ; 0 

对 图 像 的 光滑 作用 增强 ， 之 后 获取 的 各 级 反射 分 量 所 表征 的 ”用 来 调整 三 个 颜色 通道 之 间 的 比重 ; C 为 各 颜色 通道 的 颜色 

细节 信息 的 尺度 越 来 越 大 。 调整 因子 ， 通 过 与 待 色 彩 恢复 的 图 像 进行 卷 积 操作 得 到 最 终 
在 图 像 增强 的 过 程 中 ， 根 据 具 体 的 图 像 增强 场景 ， 其 对 ”的 输出 图 像 , 在 本 文 实验 中 ,颜色 恢复 相关 的 参数 设置 如 下 : 

图 像 不 同 尺度 细节 增强 的 要 求 不 尽 相 同 。 如 对 于 低 照 度 图 像 B=1.5, y-65, 0-02. 

增强 的 场景 ， 在 设计 增强 参数 时 应 适当 增加 图 像 表 征 图 像 大 针对 多 尺度 Retinex 算法 增强 后 的 图 像 存在 的 色彩 失真 

尺度 细节 的 反射 分 量 的 增强 参数 ， 对 于 高 动态 范围 图 像 的 增 问题 ， 文 献 [19] 提 出 了 色彩 恢复 的 方法 。 该 方法 的 核心 是 保 


强 ， 在 设计 时 应 适当 增加 中 尺度 和 小 尺度 细节 反射 分 量 的 增 ” 持 图 像 三 通道 颜色 之 间 的 比例 ， 并 通过 颜色 调整 因子 实现 色 
强 参 数 ， 以 使 图 像 的 纹理 细节 得 到 有 效 恢复 ， 对 于 细节 增强 ” 彩 的 恢复 和 增强 。 本 文 在 此 选用 的 公式 和 相关 参数 设置 参照 
(details exaggeration.) 的 任务 ,应 主要 通过 控制 图 像 小 尺度 文献 进行 设置 。 事 实 上 ， 在 实验 的 过 程 中 ， 通 过 控制 变量 法 
细节 反射 分 量 的 增强 参数 。 同 时 ， 考 虑 图 像 中 噪声 点 等 小 尺 0 对 三 个 参数 进行 分 析 ， 例 如 对 于 7， 让 其 在 [4,10] 间 进行 以 
度 信息 会 反映 在 一 级 反射 分 量 上 ， 为 了 避免 噪声 的 方法 ， 通 。 0.1 为 步 长 进行 调整 ， 结 果 表 明 ， 对 于 本 文 算法 及 低 照度 场 
常 将 设置 为 1; qf os Hos 的 取 值 范围 设置 为 [1, 3]， 在 处 理 。 景 ， 文 献 [19] 给 出 的 参数 设置 能 够 在 大 多 数 条 件 下 取得 比较 
低 照 度 图 像 时 ， 本 文 将 其 分 别 设置 为 1.5 和 1.3。 在 对 反射 分 ”好 的 效果 ， 表 现在 图 像 的 整体 照度 更 加 均衡 ， 图 像 信息 灶 也 
量 进行 增强 时 ， 通 常 在 对 数 域 中 将 三 个 反射 分 量 同时 进行 处 。”” 越 高 。 
理 ， 即 3.4 算法 的 形式 化 描述 
f (F,,F,, F,) = exp[x, *In(F) Q2) 本 文 算法 的 流程 为 : 

x, * In (F,) + x * In (F,)] a) 按照 式 O) — (11) 将 图 像 从 RGB 颜色 空间 转换 到 


dt 


对 照度 分 量 的 增强 ， 应 在 提升 图 像 整体 亮度 的 同时 对 高 。 HSI 空间 ， 并 提取 亮度 分 量 7 进行 处 理 ; 
亮度 区 域 进行 抑制 ， 为 此 本 文采 用 S 型 曲线 函数 对 图 像 的 三 b) 根据 本 文 在 2.2 节 提 出 基于 双边 滤波 的 图 像 多 级 分 解 
级 照度 分 量 进行 处 理 ， 即 算法 及 式 (7)， 对 亮度 通道 7 进行 分 解 ， 获 取 反 映 图 像 大 尺 
HE. RENI 度 信 息 j 级 反射 4 量 FS, 最 终 JHE 
E b EAE, 23) 度 、 中 尺度 和 小 尺度 信息 的 各 级 反射 分 及 最 终 的 


E EL; 
其 中 : a 是 照度 调整 系数 ， 其 取 值 范 围 为 [0.1, 2.0]; < 是 防 c) 为 增强 各 级 反射 分 量 的 细节 表达 能 力 ， 按 照 式 (21) 
止 除数 为 零 的 微小 扰动 变量 。 图 6 展示 了 不 同 的 4 对 应 的 函 ~ (220 对 各 级 图 像 各 级 反射 分 量 进 行 增 强 ， 获 得 增强 后 的 
数 图 像 。 低 照度 图 像 增强 要 求 对 图 像 大 部 分 区 域 进行 增强 ， 反射 分 量 f PFF) ; 
而 对 少 部 分 高 亮 区 域 进行 抑制 ， 因 此 建议 将 4 设置 为 0.5。 d) 为 增强 低 照度 图 像 的 整体 照度 分 布 ， 按照 式 (23) 对 
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3 级 照度 分 量 进行 增强 ， 获 取 增 强 后 的 照度 分 量 L) ; 相关 算法 。 算 法 涉及 的 参数 使 用 相关 论文 推荐 的 进行 设置 。 
e) 按照 公式 (1) CD 将 增强 后 的 照度 分 量 和 反射 分 量 实验 中 所 有 算法 均 基 于 Matlab 软件 平台 实现 , 同时 为 提升 双 

进行 融合 ， 获 得 增强 后 的 亮度 分 量 六 ; 边 滤波 的 效率 , 使 用 Guarnieri 等 人 (5 提出 的 加 速 的 双边 滤波 


f) 按照 式 (12) ~ 20) 将 图 像 从 HSI 颜色 空间 转换 下 方法 。 
RGB 颜色 空间 ， 并 根据 式 (250 (26) 对 增强 后 的 图 像 进行 ”4.1 主观 评价 
颜色 调整 ， 获 得 最 终 增 强 图 像 。 为 了 衡量 本 文 算法 的 有 效 性 ， 选 择 经 典 的 低 照 度 图 像 进 

行 处 理 。 图 7 展示 了 在 阴 天 环境 下 获取 的 野外 照片 ， 从 图 中 
4 。 实验 结果 和 分 析 可 以 看 出 图 像 整体 两 地 比较 低 , 左 侧 部 分 很 难 辨认 图 像 细节 。 

为 了 验证 本 文 所 提 算 法 的 有 效 性 ， 选 择 了 多 组 低 照 度 图 经 过 算法 处 理 , 都 有 效 的 提升 了 图 像 的 整体 视觉 效果 。 其 中 ， 
像 进行 实验 对 比 。 选 择 具 有 代表 性 的 六 种 低 照 度 图 像 增强 算 CLAHE 算法 有 效 的 提升 了 图 像 的 对 比 度 ， 但 同时 放大 了 图 
法 进行 对 比 ， 包 括 Zuiderveld 等 人 Bl 提出 的 受 限 对 比 度 自 适 像 中 右上 方 天 空 区域 的 噪声 并 降低 了 亮度 BPDHE 算法 对 
应 直方 图 均衡 化 算法 (contrast limited adaptive histogram 图 像 的 增强 效果 有 效 , 左 侧 部 分 仍然 难以 辨识 。 Dong 算法 是 
equalization, CLAHE), Ibrahim 等 人 05 提 出 的 亮度 保持 动态 基于 伪 雾 图 实现 的 低 照 度 增强 算法 ， 算 法 对 图 像 整体 的 处 理 
直方 图 均衡 化 算法 (brightness preserving dynamic histogram 比较 好 ,但 在 图 像 的 边缘 区 域 过 渡 不 自然 ,影响 了 整体 效果 。 
equalization, BPDHE), Dong 等 人 09 提 出 的 基于 伪 雾 图 清晰 MSRCR 算法 、ENR 算法 和 NPEA 算法 都 属于 Retinex 相关 
化 的 算法 , NASA 中 心 提出 的 带 有 颜色 恢复 的 多 尺度 Retinex 算法， 但 MSRCR 算法 处 理 后 的 图 像 存在 较为 明显 的 色彩 偏 
算法 CMSRCRO, Shin 等 人 四 提出 的 ENR 算法 ，Wang 等 人 —— 22; ENR 算法 对 图 像 的 细节 处 理 较 好 ， 但 对 于 低 照 度 图 像 存 
07 提 出 的 NPEA 算法 。 其 中 ，CLAHE 算法 和 BPDHE 算法 在 增强 不 足 的 问题 ，NPEA 算法 处 理 后 的 图 像 照 度 整体 比较 
属于 对 比 度 增 强 相 关 算 法 , Dong 的 算法 属于 基于 伪 雾 图 清晰 均衡 ， 但 图 像 整体 仍然 偏 暗 。 而 本 文 算 法 能 够 在 提升 图 像 整 
化 的 算法 ,MSRCR 算法 .ENR 算法 和 NPEA 算法 属于 Retinex ” 体 照 度 的 同时 有 效 保持 图 像 的 纹理 细节 ， 同 时 色彩 丰富 。 


(a) 原 图 (c) BPDHE 


(e)MSRCR (gNPEA (本 文 算法 


图 7 低 照度 图 像 处 理 算 法 处 理 效果 对 比 图 


Fig.7 The comparison of low-illumination image processing algorithms 


— (eMSRCR (DENR (g)NPEA 


(h) 本 文 算法 


图 8 低 照 度 图 像 处 理 算法 处 理 效果 对 比 图 
Fig.8 The comparison of low-illumination image processing algorithms 

BI 8 3X6 Ej NASA 兰 利 研究 中 心 (Langley Research Center) 增强 图 像 照度 , 甚至 降低 了 图 像 的 整体 亮度 ; Dong 算法 在 图 

官网 , 是 常用 的 图 像 增强 对 比 素材 ,由 图 中 可 以 看 出 , CLAHE ” 中 白 塔 的 边缘 存在 明显 的 黑 边 :; MSRCR 算法 处 理 后 图 像 整 

算法 处 理 后 图 像 照度 分 布 更 加 无 序 ，BPDHE 算法 未 能 有 效 ARE, 即 存在 色 偏 ENR 算法 和 NPEA 算法 对 图 像 的 增 
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强 存在 明显 不 足 。 坟 
42 客观 评价 

1) 所 给 方法 的 评价 指标 

为 客观 评价 算法 结果 ,选取 局 部 标准 差 平 
Local Standard Deviation，ALSD )， 信 "m ( 
平均 梯度 (Average Gradient, AG) 作为 评价 
ALSD 首先 将 图 像 进行 


王 oW 


等 : 多 级 分 解 的 Retinex 低 照 度 图 像 增 强 算 法 


取 分 块 图 像 的 标 ; 


ALSD 计算 公式 为 


- I; 
25x25 4 


X OD 中 M 表示 将 原 图 像 分 成 的 块 数 ， 


“均值 (Average 
Entropy)， 以 及 
指标 。 

分 块 (本 文中 分 块 尺寸 为 25x25)， 
住 差 平 均值 作为 图 像 对 比 度 评 价 指标 。 


Q7) 


Q8) 


有 表示 第 i 个 


分 块 上 第 j 个 位 置 上 的 像素 值 ， 表示 第 i 个 分 块 的 平均 像 
素 值 。 
zi E48 FI -T- HE Set PUR ILI SE CREE, dei ERRAK HH 
信息 越 丰富 ， 其 定义 为 
Entropy — X px)lg p). (29) 
X (29) 中 Pa) 表示 像素 值 为 的 概率 ，M 表示 图 像 像 
素 的 总 个 数 。 
平均 梯度 ， 反 映 图 像 细 节 对 比 的 表达 能 力 ， 它 是 图 像 清 
晰 度 的 重要 表征 。 其 计算 公式 为 
AG-A LEY, (30) 
X (30) 中 ，V; 表示 沿 i 方 向 (x 表示 水 平方 向 、}y 表 
示 垂 直方 向 ) 上 的 梯度 图 像 在 位 置 J 处 的 梯度 值 。 
表 1 展示 了 使 用 局 部 标准 差 平 均值 、 离 散 炉 及 平均 梯度 


对 低 照度 图 像 增强 算法 的 对 上 


离散 粹 和 局 部 标准 差 平均 值 两 项 


表明 本 文 算法 能 够 有 效 
晰 度 。 对 于 平均 梯度 这 一 指标 
也 算法 仍然 比较 


但 相 比 于 其 1 


这 与 


结果 。 对 于 图 


兽 强 低 


7， 本 文 算法 在 
都 取得 了 比较 高 的 得 分 ， 


照度 图 像 的 信息 


丰富 程度 和 清 


, 本 文 算法 不 如 MSRCR 算法 ， 


高 ， 但 从 主观 评 
算法 的 结果 存在 明显 的 色彩 偏差 和 照度 不 均衡 。 
ik, CLAHE 算法 和 NPEA 算法 的 主观 评价 
主观 评价 的 结果 相 一 致 。 


价 中 ，MSRCR 


除了 本 文 算 


结果 紧 随 其 后 
事实 上 ，CLAHE $i 


通过 调整 图 像 的 直方 图 分 布 实现 图 像 增强 ， 世 


能 有 效 改 


REE E, 
H (b) 


所 出 现 的 天 空 


但 容易 造成 图 像 存在 严重 的 过 


NPEA 算法 注重 图 像 自然 色 


。 对 于 图 


区 域 的 过 暗 ， 


兽 强 ， 例 如 
以 及 边缘 处 的 光 坚 效 应 。 


, 彩 的 恢复 ， 但 对 图 


看 取得 较 高 的 得 分 。 


出 的 基于 伪 雾 


| 


算法 


8， 本 文 算法 同样 在 离散 粹 


和 局 部 标准 差 


同时 相对 梯度 方 
. CLAHE 算法 和 MSRCR 算法 。 这 


面 不 如 Dong 


p 
Eal 
医 


主观 评价 的 乡 


5 ”结束 语 
针对 低 照 度 图 像 增强 算法 存在 的 问题 ， 提 


于 本 文 算 法 更 加 
像 能 够 使 得 表征 图 像 信息 丰 
高 。 综 上 所 述 ， 本 文 算法 能 够 
和 清晰 度 ， 同 时 对 图 像 对 比 度 # 


侧重 于 不 同 尺度 细节 信息 的 增强 ， 


富 程度 的 离散 焙 指 标 比 较 


吉 果 相 一 致 。 


效 提升 图 像 的 信息 丰富 程度 
的 提升 也 较为 合理 ， 


与 4.1 节 


出 了 多 级 分 解 


的 Retinex 低 照度 图 像 增 强 模 算法 。 本 文 的 主要 结论 包括 : 


a) 通过 对 现 有 双边 


滤波 性 能 的 分 析 , 提出 基于 双边 滤波 


的 Retinex 多 级 分 解 ， 从 而 获得 表征 不 同 尺度 信息 的 多 级 有 反 


射 分 量 和 最 终 的 照度 分 量 。 
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b ) 通 过 使 用 带 有 不 同系 数 的 指数 函数 对 分 解 得 到 的 多 级 
反射 分 量 进行 增强 ， 能 够 有 效 控制 图 像 各 尺度 信息 的 表达 能 
力 ， 同 时 有 效 增 强 图 像 的 边缘 细节 等 高 频 信 息 。 通 过 使 用 S 
型 曲线 函数 对 分 解 得 到 的 最 重 的 照度 分 量 进行 处 理 ， 能 够 在 
是 升 图 像 整 体 亮 度 的 同时 抑制 高 亮 ~ 最 后 ， 为 了 减少 
增强 图 像 的 色彩 偏差 和 失真 ， 使 用 颜色 恢复 函数 对 结果 进行 
后 处 理 ， 得 到 最 终 输出 图 像 。 

c) 实验 结果 表明 相 比 于 现 有 的 低 照度 图 像 ， 本 文 算法 
能 够 更 加 有 效 的 提升 图 像 的 清晰 度 `, 对比度 和 信息 丰富 
展现 了 本 文 算法 的 优越 性 。 


T 


程度 ， 


dl 低 照度 图 像 增 强 算法 客观 评价 指标 比较 
Table 1 Comparison of objective evaluation for low-illumination 
image processing algorithms 
示例 指标 原 图 像 


CLAHE BPDHE Dong MSRCR ENR NPEA 本 文 算法 


DE 7.0869 7.1806 6.7536 60.8796 7.4921 7.5885 7.1043 7.7185 


Fig.7 ALSD 15.4116 23.504 19.2311 28.2269 25.254 18.8145 22.654 36.7398 
平均 梯度 5.4102 10.9322 7.8767 11.7013 12.7391 6.9651 7.8339 12.0331 


DE 6.8651 7.5107 7.2926 7.2733 7.5316 7.6441 7.1657 7.5873 


Fig.8 ALSD 16.8432 33.8842 21.8602 35.5895 34.4887 28.3034 28.7496 44.791 
平均 梯度 5.8608 13.4004 7.8712 15.5982 13.9417 10.2034 10.48 12.7194 
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